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Using 449 million BB pairs collected with the Belle detector at the KEKB asymmetric-energy 
e + e~ collider, we observe clear signals for B + — * K°K + and B° — > K°K° decays with 5.3 a and 6.0 a 
significance, respectively. We measure the branching fractions B(B + — > K°K + ) — (1.22lQ'28io'i6) x 
10" 6 and B(B° -> K°K°) = (0.87t°'il ± 0.09) x 10~ 6 , and partial-rate asymmetries A C p(B+ -> 
K°K+) = 0.13±g;^±0.02 and A C p(B° -> ~K a K°) = -0.58^;^ ±0.04. From a simultaneous fit we 
also obtain B(B + -» K°n + ) = (22.8±g;f ± 1.3) x 10" 6 and A C p(B + -> K°tv + ) = 0.03 ± 0.03 ± 0.01. 
The first and second error in the branching fractions and the partial-rate asymmetries are statistical 
and systematic, respectively. No signal is observed for B — > K + K~ decays, and for this branching 
fraction we set an upper limit of 4.1 x 10~ 7 at the 90% confidence level. 

PACS numbers: 13.25.Hw, 11.30.Er, 12.15.Hh, 14.40.Nd 



All B — > Kir, irir decays have now been observed [J, 
0, an d direct CP violation has been established in 

The measurements of these hadronic 
b — ► s and b — > u transitions have provided essential in- 
formation for our understanding of B decay mechanisms, 
and are probes for possible new physics. What remains 
are the B — > KK modes, for which measurements with 
good accuracy are needed. Some measurements for these 
modes were reported by the Belle and BaBar collabora- 
tions 0,0. _ 

In this Letter we report the observation of B° — > K°K° 
and B+ -> K°K+ decays @. These decays are ex- 
pected to be dominated by the loop-induced b — > dss 
process (called a b — > d penguins) shown in Fig. [Ha). 
When compared with the b — > sdd penguin-dominated 
B° — > K°tt + decay, these modes are expected to be sup- 
pressed by a factor of roughly 1/20, with branching frac- 
tions expected at the 10~ 6 level fiot fTl| . We also search 
for B° — > K + K~ , which, at lowest order, arises from a 
b — » u PF-exchange process (Fig. [ljb)) or from final-state 
interactions [12j. 

For the decay modes with significant signal, we also 
measure the partial-rate asymmetry, 



A, 



CP 



N(B -» /) - N(B -> /) 



(1) 



N(B^f) + N(B^f) 
where / denotes K°K°, K°K+ or K n + . Direct CP 



violation is expected to be sizable in B° — > K a K° 
and B + — > K a K + decays [loj . while mixing-dependent 
CP violation can be measured in B° — > K°K a (and 
K+K~) [ll|. 



This analysis is based on a sample of (449.3 ± 5.7) 
xlO 6 BB pairs collected with the Belle detector at the 
KEKB e + e~ asymmetric-energy (3.5 on 8 GeV) col- 
lider The production rates of B + B~ and B°B° 
pairs are assumed to be equal. Throughout this paper, 
the inclusion of the charge-conjugate decay is implied, 
unless explicitly stated otherwise. 

The Belle detector is a large-solid-angle magnetic spec- 
trometer that consists of a silicon vertex detector (SVD), 
a 50-layer central drift chamber (CDC), an array of aero- 
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FIG. 1: The b -» d penguin diagram (a) for B + -> K°K + 
and B° — > K°K° modes, and b — » u W-exchange diagram (b) 
for B° -> K+K ~ decay. 
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gel threshold Cherenkov counters (ACC) , a barrel-like ar- 
rangement of time-of- flight scintillation counters (TOF), 
and an electromagnetic calorimeter comprised of CsI(Tl) 
crystals located inside a superconducting solenoid coil 
that provides a 1.5 T magnetic field. An iron flux-return 
located outside the coil is instrumented to detect K\ 
mesons and to identify muons. The detector is described 
in detail elsewhere [14j . Two different inner detector con- 
figurations were used. For the first sample of 152 million 
BB pairs (Set I), a 2.0 cm radius beampipe and a 3- 
layer silicon vertex detector were used; for the latter 297 
million BB pairs (Set II), a 1.5 cm radius beampipe, a 4- 
layer silicon detector and a small-cell inner drift chamber 
were used 



15] 



Primary charged tracks are required to have a distance 
of closest approach to the interaction point (IP) of less 
than 4 cm in the beam direction (z) and less than 0.1 
cm in the transverse plane. Charged kaons and pions are 
identified using dE/dx information from the CDC and 
Cherenkov light yields in the ACC, which are combined 
to form a K-n likelihood ratio, 1Z(K/tt) = Ck/{Ck + 
Lit), where Ck is the likelihood that the track is a 
kaon (pion). Charged tracks with 1Z(K/tt) > 0.6 (< 0.4) 
are regarded as kaons (pions) for B + — » K°K + (B + — > 
K°ir + ) decays. A tighter requirement, 1Z(K/tt) > 0.9, 
is used for the B° — ► K + K~ selection due to the large 
background from B — > K + ?r~. Furthermore, in all de- 
cay modes we reject charged tracks consistent with an 
electron hypothesis. 

Candidate K° mesons are observed as K s — > tt + tt~ 
decay with the branching fraction taken from Ref. [lij ]. 
We pair oppositely-charged tracks assuming the pion hy- 
pothesis and require the invariant mass of the pair to be 
within ±18 MeV/c 2 of the nominal Kg mass. The inter- 
section point of the 7r + 7r~ pair must be displaced from 
the IP. 

Two variables are used to identify B candi- 
dates: the beam-energy constrained mass, M^ c = 



P*b/c 2 , and the energy difference, AE 



Et 



^boam: where £ bcam 



is the run-dependent beam 
energy and Eg and p* B are the reconstructed energy 
and momentum of the B candidates in the center-of- 
mass (CM) frame, respectively. Events with Mhc > 5.20 
GeV/c 2 and \AE\ < 0.3 GeV are selected for the analy- 
sis. 

The dominant background is from e + e~ —> qq {q = 
u, d, s, c) continuum events. We use event topology to 
distinguish the spherically distributed BB events from 
jet-like continuum background. We combine a set of 



modified Fox- Wolfram moments [17| into a Fisher dis- 
criminant. Signal and back grou nd likelihoods are formed, 
based on a GEANT-based [HI Monte Carlo (MC) simu- 
lation, from the product of the probability density func- 
tion (PDF) for the Fisher discriminant and that for 
the cosine of the angle between the B flight direction 
and the positron beam. Suppression of continuum is 



achieved by applying a requirement on a likelihood ra- 
tio TZ = £ S ig/(£ si g + £ q q), where £ sig {C q q) is the signal 
(qq) likelihood. Additional continuum background sup- 
pression is achieved through use of a £?-flavor tagging 
algorithm [19], which provides a discrete variable indi- 
cating the flavor of the tagging B meson and a quality 
parameter r, with continuous values ranging from (for 
no flavor-tagging information) to 1 (for unambiguous fla- 
vor assignment). Events with a high value of r are con- 
sidered well-tagged and hence arc unlikely to have orig- 
inated from continuum processes. We classify events as 
well-tagged (r > 0.5) and poorly-tagged (r < 0.5), and 
for each category of Set I and Set II we determine a con- 
tinuum suppression requirement for TZ that maximizes 

Here, denotes 



N c - P 
qq ' 



the value of Ng/^Ng 

the expected signal yields based on MC simulation and 
the branching fractions of the previous measurement 0] , 
and JV^ P denotes the expected qq yields as estimated 
from sideband data (M bc < 5.26 GeV/c 2 and \AE\ < 0.3 
GeV). _ 

Background contributions from other T(4S) — > BB 
events are investigated with a large MC sample that in- 
cludes events from b — ► c transitions and charmless B 
decays. After all the selection requirements, no BB back- 
ground is found for the B° — > K°K° mode. Due to K-tt 
misidentification, large B° — > K + tt~ and B + — > K°tt + 
feed-across backgrounds appear in the B° — ► K + K~ and 
B + — ► K K + samples, respectively. At low AE values, a 
small background contribution from charmless B decays 
- mainly from B+ -> K+7r°K (B° -> K*-n+) — is 
found for the B+ -> K°K+ {B+ -> iv°7r+]_decay mode. 
However, experimentally the B + — > K + ir°K° mode only 
has an upper limit. 

The signal yields are extracted by performing extended 
unbinned two-dimensional maximum likelihood (ML) fits 
to the (Mbc, AE) distributions. The likelihood for each 
mode is defined as 

C = eM-^Nhk^U^^Hk^ ( 2 ) 

l,k,j i l,k,j 

Kkj = \[l~q l -Acp^P^iM^AE*), (3) 

where i is the event identifier, I indicates Set I or Set II, 
k distinguishes the two r categories and j runs over all 
components included in the fitting function — one for 
the signal and the others for continuum, feed-across and 
charmless B backgrounds. N^kj represent the number 
of events, and Pi,k,j(M£ c , AE 1 ) are the two-dimensional 
PDFs, which are the same in the two r categories for 
all fit components except for the continuum background. 
The parameter q indicates the B-meson flavor: q = 
+1 (-1) for B+ and B° (B~ and B°). Unlike K°K+, 
the K°K° and K + K~ channels are not self-tagged and 
the B meson flavor must be determined from the other 
B. To account for the effect of BB mixing and imper- 
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feet tagging, the term Acp for the signal in Eq. (3) has 
to be replaced by A C p(l - 2xd)(l - 2wk), where Xd = 
0.188 ± 0.003 [16| is the time-integrated mixing parame- 
ter. The Xd value of continuum events is set to zero. The 
wrong-tag fraction Wk, which depends on the value of r, 
is determined from a high statistics sample of self-tagged 
B° -> D*-n+,D*-p+ and D^~l+Ui events [H. 

All the signal PDFs (Pi,k,j=signa.l(M hc , AE)) are pa- 
rameterized by smoothed two-dimensional histograms 
obtained from correctly reconstructed signal MC based 
on the Set I and Set II detector configurations. Signal 
MC events are generated with the PHOTOS [2(| sim- 
ulation package to take into account final-state radia- 
tion. Since the Mb c signal distribution is dominated by 
the beam-energy spread, we apply small corrections to 
the signal peak position and resolution determined using 
B+ -> ~D q tt+ from data (D° -> K^tt+tt- is used for the 
K°K° mode, while D° — > K + ir~ is used for the other 
three modes) with small mode-dependent corrections ob- 
tained from MC. The resolution for the AE distribu- 
tion is calibrated using the invariant mass distributions 
of high momentum (pLab > 3 GeV/c) D mesons. The 
decay mode D° — > K + n~ is used for B° — > K + K~, 
D+ -> K° s ir+ for B+ -> K°ir+ and U° -> ivgrr+TT- for 

The continuum background PDF is described by a 
product of a linear function for AE and an ARGUS 
function, f(x) 



M hc c 2 /E } 



cam 



xyl — x 2 exp[— £(1 — x 2 )}, where x 
The overall normalization, AE 
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slope and ARGUS parameter £ are free parameters in 
the fit. The background PDFs for charmless B de- 
cays for the K°tt + and K°K + modes are both mod- 
eled by smoothed two-dimensional histograms obtained 
from a large MC sample. We also use smoothed two- 
dimensional histograms to describe the feed-across back- 
grounds for the K + K~ (K°K + ) mode, since the back- 
ground K + ir~ (K°ir + ) events have (Mb c ,AE) shapes 
similar to the signal, except for the AE peak positions 
shifted by ~ 45 MeV. We perform a simultaneous fit for 
B + — ► K°K + and B + — > i4T°7r + , since these two decay 
modes feed into each other. Because the branching frac- 
tion of B° — > K + K~ is small, we also treat the yields of 
K + ir~ feed-across events as free parameters in the fit. 

When likelihood fits are performed, all the Ni^j are 
allowed to float except for the feed-across backgrounds 
in the K°K + and K°n + modes. The M hc and AE pro- 
jections of the fits are shown in Fig. O The branch- 
ing fraction in each mode is calculated by dividing the 
efficiency-corrected total signal yield by the number of 
BB pairs. In Table [H a sum of fitted signal yields and 
the average efficiency are listed. 

The fitting systematic errors include the signal PDF 
modeling, which we estimate from the deviations after 
varying each parameter of the signal PDFs by one stan- 
dard deviation in the calibration factors, and the model- 
ing of the charmless B background. Since the AE val- 



ues of the charmless B events are typically less than 
—0.12 GeV, the systematic error due to the modeling 
of the charmless B background is evaluated by requir- 
ing AE > —0.12 GeV. At each step, the yield deviation 
is added in quadrature to provide the fitting systematic 
errors, and the statistical significance is computed by tak- 
ing the square root of the difference between the value of 
— 2 ln£ for the best fit and that for zero signal yield. For 
B + — > K°K + , systematic uncertainty is included in the 
significance calculation by varying the feed-across back- 
ground (which is the dominant uncertainty) by la in the 
direction that lowers the significance. For the other de- 
cay modes, the effect of systematic uncertainty on the 
significance is negligible. 
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FIG. 2: M bc (left) and AE (right) distributions for B -> 
K+K~, K°K+, K°ty+ and K°K° candidates. The his- 
tograms show the data, while the curves represent the var- 
ious components from the fit: signal (dot-dashed), continuum 
(dashed), charmless B decays (hatched), feed-across back- 
ground from misidentification (dotted), and sum of all com- 
ponents (solid). The Mb c and AE projections of the fits 
are for events that have \AE\ < 0.06 GeV (left) and 5.271 
GeV/c 2 < Mbc < 5.289 GeV/c 2 (right), respectively. 



The MC-data efficiency difference due to the require- 
ment on the likelihood ratio 1Z is investigated using the 
B+ -► D°ir+ (D° -► K^tt+it- for K°K° and D a -> 
K + ir~ for the others) samples. The systematic error 
due to the charged-track reconstruction efficiency is esti- 
mated to be 1% per track using partially reconstructed 
D* events. The systematic error due to H(K/ir) selection 
is 1.3% for pions and 1.5% for kaons, respectively. Due 
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TABLE I: Fitted signal yields, product of efficiencies and sub-decay branching fractions (B s ), branching fractions, significance 
(£), and partial-rate asymmetries for individual modes. The first and second error in the branching fractions and the partial-rate 
asymmetries are statistical and systematic, respectively. 



Mode 


Yield 


Eff.x6 s (%) 


£(HT fa ) 


E(ct) 


A C p 


K+K- 


o c+5.0 
z -°-3.7 


6.18 


Q.09ig;i| ±0.01 (< 0.41) 


0.6 




K°K+ 


36.6+g 


6.72 


i oo +0.32+0.13 
1 -"-0.28-0. 16 


5.3 


0.13la24 ±0.02 


K°TV + 


1 252 


12.21 


22.8lo. 7 ± 13 


53.1 


0.03 ±0.03 ±0.01 


K°K° 


23.0+1!; 


5.89 


0.87+0 20 ± 0.09 


6.0 


-0.58+Q gg ± °- 04 



TABLE II: Summary of systematic errors, given in percent. 





K+K- 


K U K+ 


A'° 71 + 


A U A U 


Signal PDF 


+ 1.3 
-1.4 


±0.2 


±0.2 


+0.5 

-0.6 


Charmless B background 


0.0 


-0.9 


-0.1 


0.0 


1Z requirement 


±0.8 


±1.4 


±1.1 


±3.3 


Tracking 


±2.0 


±1.0 


±1.0 


0.0 


lZ(K/ir) requirement 


±3.9 


±1.5 


±1.3 


0.0 


Kg reconstruction 


0.0 


±4.9 


±4.9 


±9.8 


# of feed-across 


0.0 


+9.4 
-11.9 


+0.2 
-0.4 


0.0 


# of BB 


±1.3 


±1.3 


±1.3 


±1.3 


Signal MC statistics 


±1.1 


±1.1 


±0.8 


±0.8 


Sum 


+4.9 
-5.0 


+ 11.0 

-13.2 


±5.5 


±10.5 



to the tighter 1Z(K/tt) selection of kaons in the K + K~ 
mode, the assigned systematic uncertainty is 1.9% per 
kaon track. The Kg reconstruction is verified by com- 
paring the ratio of D + — > Kgn + and D + — » K~tt + tt + 
yields with the MC expectation. We vary the yields of 
feed-across background by ±1<7 to check the effect from 
the constraint on the feed-across background. Possible 
systematic uncertainties due to the description of final- 
state radiation have been studied by comparing the lat- 
est theoretical calculations with the PHOTOS MC Q. 
These uncertainties were found to be negligible and thus 
no systematic error is assigned due to PHOTOS. The 
systematic error due to the uncertainty in the total num- 
ber of BB pairs is 1.3% and the error due to signal MC 
statistics is in the range 0.8 - 1.1%. The final systematic 
errors are obtained by quadratically summing the errors 
due to the reconstruction efficiency and the fitting sys- 
tematics. The summary of the systematic errors is shown 
in Table HU 

The detector bias is the dominant systematic error 



for A CP (B+ -► K°K+) and A C p(B^ 



"); the 



systematic uncertainties evaluated from the partial rate 
asymetry of the continuum background are 0.02 and 0.01 
for these two modes, respectively. The systematic errors 
for A C p(B" -> K°K°) are estimated by varying the fit- 
ting parameters by ±lcr. We include also the errors due 
to Wk , Xd and tag-side interference [23| and obtain a total 
systematic error of 0.04. 

In summary, using a data sample with 449 mil- 



lion BB pairs, we observe B + — » K°K + and B° — » 
~K°K° with branching fractions B{B+ -> ~K a K+) = 
(1.22+°;^+° 1 1 3) x 10- 6 and B(B° - K°K°) = 
(0.87to;2o±0.09) x 10~ 6 . The corresponding partial-rate 
asymmetries are A CP {B+ -> K°K+) = 0.13±g;l| ± 0.02 
and A C p{B° -> K°K°) = -0.58 ± 0.04. In ad- 
dition, we improve the measurements of the branching 
fraction and partial-rate asymmetry for the decay B + — > 
Ar°7T+: B{B+ -> A'°7r+) = (22.8±g;f ± 1.3) x 10~ 6 and 
A C p(B+ -> K°tt+) = 0.03 ± 0.03 ± 0.01. Our measure- 
ments are consistent with previous results [H, 0, S] . The 
new results, except for B(B° — > K + K~) and A C p(B — > 
AT°A' ), have better precision than previously measured 
values. Our result s ag ree with some theoretical predic- 
tions U , 24 , l25l I2II ] . No signal is observed for 
B° — > A'+A" - , and we set an upper limit of 4.1 x 10~ 7 
at the 90% confidence level using the Feldman-Cousins 
approach [27]. 
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